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ABSTRACT 



Aims. The existence of an anti-correlation between the Equivalent Width (EW) of the neutral narrow core of the iron Ka emission 
line and the 2-10 keV luminosity (the so-called 'X-ray Baldwin' or Twasawa-Taniguchi' effect) has been debated in the last years. 
We aim at testing this claim on the largest catalogue of radio quiet AGN high-quality X-ray spectra ever published. 
Methods. The final sample comprises 157 objects. We search for a relation of the iron line EW not only with the X-ray luminosity, 
but also with the Black Hole mass, the Eddington ratio and the cosmological distance. The data presented here were analyzed homo- 
geneously, all spectra are from the same instrument and with high Signal-to-Noise Ratio. 

Results. A linear censored fit on the EW versus 2-10 keV luminosity is highly significant and yields log(EW Fe ) = (1.73 ± 0.03) + 
(-0. 17 ± 0.03) log(Lx,44), where EWfe is the EW of the neutral iron Ka line in eV and L x ^ is the 2-10 keV X-ray luminosity in units 
of 10 44 erg s -1 . The anti-correlation with the Eddington ratio is also very significant, while no dependence of the iron EW on the BH 
mass is apparent. 
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1. Introduction 



iBaldwinl i 1977b first reported a significant anti-correlation be- 
tween the EW of the [Civ] A 1549 line and the UV luminos- 
ity in quasars. The so-called 'Baldwin- effect' was later on con- 
firmed on more robust basis (see e.g. iKinnev et al.lll990l) and 
similar anti-correlations with luminosity were found for other 
emission fines prod uced in the Broad Line Region (BLR) (e.g. 
iDietrich et al.l l2002). However, the physical explanation for this 
phenomenon is still unclear and several possibilities are gener- 
ally invoked in literature, including a change of the ionizing con- 
tinuum and gas metallicity with luminosity dKorista et al.l ll998) 
or a luminosity-dep endent covering factor and ionization pa- 
rameter of the BLR (iMushotzkv & Ferlandll 19841) . On the other 
hand, some studies claim that the primary physical parameter 
which drives the Baldwin effect may be the accretion rate, in- 
stead o f the luminosity. 

In 119931 llwasawa & Taniguchil presented a similar anti- 
correlation between the iron Ka emission line and the 2-10 
keV luminosity, according to Ginga observations of 37 AGN. 
This Twasawa-Taniguchi effect' (TT effect' from now on, also 
known as 'X-ray Baldwin effect') was then found in a larger 
sample of objects observed by XMM-Newton, giving a rela- 
tion between luminosity and Equivalent Width (E W) of the 
narrow core of the Fe Ka as EW oc £-0.i7±o.08 (Pag e et all 



2004 ). The significance o f this effect was questioned by 



Jimenez-Bailon et"aT] d2005l) in their analysis of XMM-Newton 



data of PG quasars, pointing out the i mportance of conta mina- 
tion from radio-loud objects. Recently. iJiang et ail (f2006) com- 
bined Chandra and XMM-Newton data and suggested that, af- 
ter excluding radio-loud objects, the anti-correlation can be at- 
tributed to variations of the continuum, while the iron line stays 



Send offprint requests to: Stefano Bianchi 
e-mail: bianchi @f is . uniroma3 . it 



constant, as ex pected if the two c omponents are produced on 
different scales. iNandra et alJ (1 19971) also reported an IT effect 
in their ASCA sample of AGN, but concluded that the effect was 
entirely due to the relativistically broadened component of the 
iron line. The presence of an anti-correlation between the iron 
relativistic line a nd the luminosity was also recently found by 
iGuainazzi et al. (2006a) on a large data set of XMM-Newton 
spectra of AGN. 

In this paper, we search for the IT effect in the largest cata- 
logue of radio-quiet AGN observed by XMM-Newton published 
so far. Radio-loud objects are excluded on the basis of the radio- 
loudness parameters R and Rx, calculated for each source from 
data in the literature. The data presented here were analyzed ho- 
mogeneously, all spectra are from the same instrument and with 
high Signal-to-Noise Ratio (SNR). 



2. The catalogue 

2. 1. X-ray data reduction and analysis 

Our catalogue consists of all radio-quiet Type 1 AGN observed 
as a main target by XMM-Ne wton, whose data a re public as of 
March 2007. Only EPIC pn (Str uder et al"] 12001 1) data were re- 
processed, with SASv7.0 (iGabriel et al.ll2004l) . For the obser- 
vations performed in Small Window mode, background spec- 
tra were generated using black-field event lists, according to 
the procedure presented in lRead & Ponmanl (120031) . In all other 
cases, background spectra were extracted from source-free re- 
gions close to the target in the observation event file. 

Source extraction radii and screening for intervals of flar- 
ing particle background were performed via an iterative process 
which lea ds to a maximiza t ion of the SNR, similarly to what de- 
scribed in iPiconcelli et al .1 d2004l) . Spectra were binned in order 
to oversample the intrinsic instrumental energy resolution by a 
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factor > 3. Moreover, spectra were further re -binned in order to 
have spectral bins with at least 25 background-subtracted counts, 
to ensure the applicability of the^ 2 statistics. 

We applied a number of criteria to filter the catalogue. First 
of all, spectra with less than 200 counts in either of the (rest- 
frame) bands of 0.5-2 and 2-10 keV were rejected, because 
they do not possess enough independent bins to be fitted with 
our models. Moreover, spectra affected by a pileup larger than 
1% were rejected. Only two objects, namely NGC 2992 and 
ESO548-G081, did not have an alternative pileup-free observa- 
tion and were therefore excluded from the catalogue. After the 
fitting procedure, all sources with a local column density (mea- 
sured in the 2-10 keV band) larger than 2 x 10 22 crrT 2 were fur- 
ther excluded from the catalogue. The exclusion of radio-loud 
objects is crucial for the analysis of this paper. Therefore, we 
collected radio data at 6 and 20 cm fr om the literature an d calcu- 
lated the radio-loudness parameter R { Stock e et alj 19921) and th e 
X-ray radio-loudness parameter Rx dTerashima & Wilsonl2 003). 
According to the standard definition, all quasars with \og{R) > 1 
were excluded, while for Seyfert galaxies, known to be on av- 
erage 'radio-louder', lo gQR) > 2.4 and log( Rx) > -2.755 were 
used as boundaries (see Panessa et al. 2 0071) . 

At the end of this selection procedure, the total catalogue 
comprises 157 radio quiet AGN, whose X-ray data are comple- 
mented by H/3 Full Width Half Maximum (FWHM) and Black 
Hole (BH) mass measures, when available in literature (the cov- 
erage is 64% and 52%, respectively). We defer the reader to 
Bianchi et al., in preparation, for the full presentation of the cat- 
alogue. When more than one observation of the same object was 
available, we used the one with the longest exposure, indepen- 
dently from its flux status, in order to avoid, as much as possible, 
any bia s for the choic e. All spectra were analyzed with Xspec 
12.2.1 dArnaudl [l996). For the aims of this paper, we restricted 
the analysis to the 2-10 keV band, adopting a baseline model 
consisting of a powerlaw plus three narrow (<x = 0) Gaussian 
lines with energies fixed at 6.4, 6.7 and 6.96 keV, corresponding 
to emission from neutral, He- and H-like iron. A relativistically 
broadened component of the line was never included in our fits. 
We will discuss its possible contribution to the EW of the narrow 
core in Sect. [4] In the following, errors correspond to the 90% 
confidence level for one interesting parameter (A^ 2 = 2.71), 
where not otherwise stated. The adopted cosmological param- 
eters are Hq = 70 km s _1 Mpc~', Ao = 0.73 and qo = (i.e. the 
default ones in Xspec 12.2.1). 



3. The IT effect 

Fig. Q] shows the IT effect in our data. As a result of the fitting 
procedure described in Appendix [Al we obtained the following 
best fit (see also Table [TJ: 

\og(EW Fe ) = (1.73 + 0.03) + (-0.17 + 0.03)log(L x ,44) (1) 

where EWf c is the EW of the neutral iron Ka line in eV and 
Lx,44 is the 2-10 keV X-ray luminosity in units of 10 44 erg s _1 . 
The Spearman's rank coefficient for this correlation is -0.33, cor- 
responding to a Null Hypothesis Probability of 4 x 10~ 5 . 

Our result is consistent with previous claim for the IT effect, 
whether radio -loud objects were inc l uded or not. In the orig- 
inal paper by Iwasawa & Taniguchi (1993), the slope of their 
anti-correlation was -0.20 + 0.03. The foll owing results ar e 
all consistent with each other: -0-17 + 0.0 8 dPage et al.ll2004l). 
-0.06 + 0.20 dJimenez-Bailon et al.1 [2005b . As for Uiang et all 



Table 1. Best fit parameters for the IT effect on the catalogue 
presented in this paper. The first line, in bold face, refers to 
the entire catalogue, while the following lines refers to Seyfert 
galaxies (S), Quasars (Q), Broad-line objects (B) and Narrow- 
line objects (N). Correlations with Eddington ratio (LD/C: 
luminosity-dependent or constant X-ray bolometric correction), 
BH mass and redshift follow. 
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Notes: Col. (l)-(2)-(3) Total number of points, measures and y upper 
limits; Col. (4)-(5)-(6)-(7) Best fit parameter and errors for the linear fit 
log(j) = A + Blog(jc); Col. (8)-(9) Mean Spearman's rank coefficient 
and Null Hypothesis probability. See AppendixfA]for details. 



(2006), they quote -0.20 + 0.04 for their whole sample and 
-0.10 + 0.05 for a sub-sample without radio-loud objects. 

We also tested if this anti-correlation still holds for sub- 
classes of our catalogue. As shown in Table [1] there is no sig- 
nificant difference between Seyfert galaxies and quasars, nor be- 
tween broad- and narrow-line objects, even if the lower statistics 
prevent us from drawing firm conclusions. 

In order to understand which is the main driver for the IT 
effect, we checked the behaviour of the iron Ka EW with re- 
spect to the Eddington ratio, the BH mass and the redshift (see 
Table[T]and Fig.|2]i. The anti-correlation with the Eddington ra- 
tio is highly significant (in particular when usin g a luminosity- 
depen dent X-ray bolometric correction, like in iMarconi et al.l 
2004, instead of a constant one). On the other hand, the cor- 
relation with z is weaker and, quite interestingly, no significant 
dependence of the iron EW on the BH mass is apparent. 

4. Discussion 

We have shown in the previous section that our data confirm 
the significance of the IT effect. Several explanations have been 
proposed in the literature for this phenomenon. 
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Iwasawa-Taniguchi effect: log(EW Fo ) = (1.73±0.03) + (-0.17+0.03) logtUjJ 




10" 10' 
2-10 keV X-ray Luminosity (10 44 , erg s' 1 ) 



Fig. 1. The TT effect': neutral iron EW against 2-10 keV X-ray luminosity of the objects in our catalogue. The anti-correlation 
between the two parameters is shown as the best fit line, whose analytical expression is reported on the top. The broken lines 
represent the combined error on the slope and normalization of the best fit. The different symbols refers to the classification of the 
objects, on the basis of their absolute magnitude and HB FWHM: NLSY, narrow-line Seyfert 1 ; BLSY, broad-line Seyfert 1 ; NCSY, 
not-classified Seyfert 1 (no HB FWHM measure available); NLQ, narrow-line quasar; BLQ, broad-line quasar; NCQ, not-classified 
quasar (no HB FWHM measure available). See text for details. 



log(EW F9 ) = (1.61±0.05) + (-0.19 + 05) log(L MI /L E0lJ ) 



log(EW Fe ) = (1 .73 + 0.04) + (-0.07 + 0.04) log(M BHjj ) 





M BH (10 M mn ) 



Fig. 2. Left panel: neutral iron EW against the Eddington ratio L[, \ / L et id, adopting a luminosity-dependent X-ray bolometric correc- 
tion. Right panel: The same as above, but against the BH mass. See caption of Fig. Q~]for details. 



The original Bal dwin effect on the JClvl /U549 line 
was interpreted by iMushotzkv & Ferlandl (Q984) in terms 
of a luminosity-dep e ndent covering factor of the BLR. 
Iwasawa & Taniguchi (1993) noted that this explanation was 
also a good one for their recently discovered X-ray Baldwin ef- 
fect, suggesting an origin of the iron Ka line in a material close 
to the BLR itself. However, even if a contribution from the BLR 
cannot be excluded for most o bjects (and in some ca ses seems to 
be the best possibility, see e.g. lBianchi et al ] l2003h. no correla- 
tion exists between the width of Fe K and the HB (Nandra 2006). 
Instead, the ubiquitous presence of an accompanying Compton 
reflection component strongly favours a sce nario where most o f 
the iron line flux comes from the torus (e.g. Bian chi et alj |2004). 

The IT effect should rather be associated to a change of the 
covering factor of a compact Compton-thick torus. Such a possi- 



bility could be related to disc-driven hydromagnetic wind mod- 
els, which predic t an increase of the op ening angle of the torus 
with luminosity (Konig i"& Karti d fl 9941) . This, in turns, would 
lead to a reduction of the covering factor of the torus, which 
could explain the observed iron line-luminosity anti-correlation. 
A side implication of this luminosity-dependent covering factor 
for the torus is a decrease of the fraction of obscured AGN with 
luminosity, which is indeed w hat recently observed ( Ued a et all 
120031 ; lLa Franca et alj 12005). It must be noted, however, that 
the latter observational result is based on Compton-thin sources 
only, while the torus is likely Compton-thick. 

Indeed, the decrease of obscured sources with luminosity can 
be explained by gravitational effects of the BH on the molec - 
ular gas in the disk of the host galaxy (Lamas tra et alj 12006). 
These authors showed that the covering factor of this large-scale, 
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log(Ratio Fe ) = (0.12±0.06) + (0.09 +0.06) logtl^ J 
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Fig. 3. The ratio between the ionized and the neutral iron lines 
against the 2-10 keV X-ray luminosity for the objects in our cat- 
alogue. See caption of Fig.[T]for details. 



Compton-thin material decreases with luminosity. Lamastra et 
al. (in prep.) investigated the iron line emission from the molec- 
ular gas in the disk, as expected from their model. As a natural 
consequence of their assumptions, the EW appears to be anti- 
correlated with the X-ray luminosity indirectly, the BH mass 
being the principal driving mechanism. However, Fig. [2] (right 
panel) seems to indicate that no anti-correlation with the BH 
mass is observed. Moreover, due to the fact that the material is 
Compton-thin, the EW calculated by Lamastra et al. (in prep.) 
can only account for less than a half of the observed iron line 
EW, at each luminosity. Therefore, it seems unlikely that their 
model can explain the observed IT effect. 

Another possibility is that the ionization state of the iron 
emitting material c hang es according to l uminosity. Even if 
iNandra et ail dl997l) and iNavakshinl {2000) proposed this sce- 
nario for the iron line component produced in the accretion 
disk, it is still possible that even the inner walls of the torus 
may become so ionized to suppress iron line emission for high- 
luminosity objects. As the ionization state rises, one would 
expect a strong contribution from He- and H-like iron (e.g. 
iMatt et aD 19961) . We tested this scenario in Fig. [3](right panel), 
where the ratio of the overall contribution from highly ionized 
iron to neutral iron is plotted against X-ray luminosity. The best 
fit between these quantities suggests indeed a weak correlation, 
but its signific ance is poo r (see a lso TableQ}. 

Recently, Jian g et all ((2006) proposed that the IT effect 
can be completely accounted for by variability effects, assum- 
ing constant iron line fluxes. This naturally produces an anti- 
correlation between X-ray luminosity and iron EW, which may 
contribute to the overall IT effect. How ever, the simulate d anti- 
correlation has a slope of -0.05 + 0.05 dJiang et al.ll2.Q06l) , much 
weaker than the observed IT effect. In any case, the amplitude 
of variability in radio-quiet AGN is generally smaller than one 
order of magnitude, while the IT effect is observed on six orders 
of magnitude in luminosity (see Fig.[T]i , making this explan ation 
unlikely, as originally noted bv llwasawa & Ta niguchi ( 1993). 

A last comment must be spent on the possibility that the 
IT effect may be completely or partly due to the relatiyistically 
broadened component of the iron line. Guainaz zi et al.l ((2006a) 
showed that a broad profile is significantly more common in low 
luminosity AGN. Their stacked spectral residuals analysis on 
=2 100 sources clearly presents a weaker red wing and narrow 



core for more luminous objects. Even though it is very difficult 
to quantitatively disentangle the contribution of each component 
to the overall EW decrease, due to the limited statistics of most 
objects, the EWs measured in our fits are likely dominated by the 
narrow core, especially when the underlying relativistic compo- 
nent is very broad. Therefore, the relativistic iron line can possi- 
bly contribute to the overall anti-correlation, but cannot account 
for most of the observed Iwasawa-Taniguchi effect, which im- 
plies an =a 80% decrease of the iron EW in four orders of mag- 
nitude in X-ray luminosity (see our best fit in Eq. (fTJ). 
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Appendix A: Linear fits on censored data 

Since we are dealing with a la rge number of upper limits for the EWs, we fol- 
lowed [Guainazzi et al. 1 2006b) to perform a linear fit, taking into account both 
upper limits and errors on the y parameter. Thei r method is an exte nsion of the 
regres sion method on censored data described by Schmitt 1 1985) and Isobe et al. 
1 1986) and can be summarized as follows. A large number of Ordinary Least 
Squares (OLS: Y vs. X) fits was performed on a set of Monte-Carlo simulated 
data derived from the experimental points according to the following rules: a) 
each detection was substituted by a random Gaussian distribution, whose mean 
is the best-fit measurement and whose standard deviation is its statistical uncer- 
tainty; b) each upper limit U was substituted by a random uniform distribution 
in the interval [0,U]. The mean of the slopes and the intercepts derived from the 
fits of each data set are our 'best fit'. Note that errors and upper limits on the x 
parameter are not considered in the procedure, although plotted in Fig.ff]|2]and 
[5] This does not affect the results significantly, given the fact that they are typi- 
cally much smaller than the errors on the EW. Moreover, the errors on the best 
fit parameters only take into account the experimental uncertainties, being the 
dispersion of these parameters within the Monte Carlo simulated data sets. Any 
intrinsic, physical scattering of the data points is not included. To test the signif- 
icance of this anti-correlation, we calculated the Spearman's rank coefficient for 
each Monte Carlo simulated data set: the mean value is taken as the Spearman's 
rank coefficient of our 'best fit' . 
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